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Quantum oscillations in graphene is discussed. The effect of interactions are addressed by Kohn's 
theorem regarding de Haas-van Alphen oscillations, which states that electron-electron interactions 
cannot affect the oscillation frequencies as long as disorder is neglected and the system is sufficiently 
screened, which should be valid for chemical potentials not very close to the Dirac point. We 
determine the positions of Landau levels in the presence of potential disorder from exact transfer 
matrix and finite size diagonalization calculations. The positions are shown to be unshifted even 
for moderate disorder; stronger disorder, can, however, lead to shifts, but this also appears minimal 
even for disorder width as large as one-half of the bare hopping matrix element on the graphene 
lattice. Shubnikov-de Haas oscillations of the conductivity are calculated analytically within a 
self-consistent Born approximation of impurity scattering. The oscillatory part of the conductivity 
follows the widely invoked Lifshitz-Kosevich form when certain mass and frequency parameters are 
properly interpreted. 



I. INTRODUCTION 

Quantum oscillations of magnetization (de Haas-van 
Alphen effect: dHvA), of conductivity (Shubnikov-de 
Haas effect: SdH) and that of Hall coefficient (Rh) 
are excellent probes for topologies of Fermi surfaces and 
properties of Fermi liquids.^ Recently interesting experi- 
ments have led to considerable insight into physical sys- 
tems su ch as graphen^ and high temperature supercon- 
ductors.'^ These oscillations fundamentally arise from 
Landau level quantization and hence their frequencies 
should be robust with respect to electron-electron inter- 
action, crystalline potential, as long as an effective con- 
tinuum theory exists, and under modest broadening of 
the Landau levels due to impurity scattering. On the 
other hand the waveform (harmonic content) or ampli- 
tudes are sensitive to many details. 

An elegant theorem for a translationally invariant 
continuum system in two dimensions and short range 
electron-electron interactions were formulated by KohiP 
to reveal the robustness of the frequencies. This is in con- 
trast to the more complex many-particle analysis of Lut- 
tinger,"^ which, in principle, can also address the wave- 
form or the amplitude of the oscillations, but appears to 
fail in two dimensions.'^ The reason behind Kohn's the- 
orem is intuitively clear. The magnetic field is never a 
small perturbation, especially in two dimensions, where 
the spectrum converts from continuous to discrete in its 
presence. Thus, the problem must be exactly diago- 
nalized with the magnetic field before considering the 
perturbative effects of short-ranged electron-electron in- 
teraction, which can in principle be well behaved, mod- 
ulo a quantum phase transition or a broken symmetry. 
Moreover, when the Landau levels are completely filled, 
the ground state is non-degenerate, resulting in a special 
stability, similar to magic nuclei. As the magnetic field 
is tuned through a sequence of filled Landau levels, the 
macroscopic state repeats itself, hence the periodicity; of 
course, this is strictly valid only if the relevant Landau 
level index, n ^ 1. 



One of the purposes of the present paper is to revisit 
Kohn's theorem in the context of graphene where SdH 
oscillation frequencies have been shown to be remark- 
ably close to what a pure noninteracting theory would 
predict.l^ We show that this is a consequence of Kohn's 
theorem despite the broken translational symmetry due 
to the crystalline potential, strengthening further the ev- 
idence of the Weyl fermionic character of the excitation 
spectra. In the process, we bring out some of the subtle 
aspects of the theorem and its connection with the more 
familiar Luttinger theorem regarding the volume of the 
Fermi surface.^' 

The second purpose is to provide an exact calcula- 
tion of the Landau levels in the presence of disorder and 
an analytical self-consistent calculation of SdH within 
self-consistent Born approximation (SCBA). A non-self 
consistent calculation was previously reporrted,!^ as well 
as a self-consistent calculation for unitary scattererd^Sl — 
perhaps such extreme strong scattering mechanism is not 
applicable to realistic graphene samples. Not only do we 
demonstrate that the oscillation frequencies are unshifted 
to a good approximation for even moderate disorder, but 
also obtain the amplitude and the waveform. The re- 
sults, when suitably interpreted in terms of a mass pa- 
rameter and a frequency scale, are similar to the widely 
used Lifshitz-Kosevich formulajf^ which, in a strict sense, 
cannot be applied in two dimensions. From our exact 
numerical calculations, we also show that for very strong 
disorder the Landau levels are shifted, but only mini- 
mally, bearing some resemblance to the unitary scatterers 
treated previously.'^ It is probably true that our results 
also imply that, to lowest order, results would remain 
unchanged when both electron-electron interaction and 
impurity scattering are considered together, at least for 
weak disorder relevant to experiments. The extension to 
longer ranged impurity scattering is straightforward and 
is not discussed. 

The plan of the paper is as follows: in Sec. II we dis- 
cuss Kohn's theorem and apply it to graphene. In Sec. Ill 
we provide exact transfer matrix and diagonalization cal- 
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culations to determine the position of the Landau levels 
and their shifts due to impurity scattering. In the second 
part of this section we calculate SdH oscillations within 
a self-consistent Born approximation (SCBA) developed 
by Andoli^lThe brief Sec. IV contains our concluding re- 
marks and there is an Appendix containing calculational 
details. 



II. KOHN'S THEOREM FOR QUANTUM 
OSCILLATIONS IN TWO DIMENSIONS 

A. Nonrelativistic fermions 

There is a precise theorem by Kohn that in a two 
dimensional continuum system electron-electron interac- 
tion cannot shift the dHvA frequency in all orders in 
perturbation theory.'^ To clarify its content, it is useful 
to reconsider it. Let the unperturbed problem be defined 
by the non-interacting Hamiltonian, Hq, in a rectangular 
box Lx X Ly. In the Landau gauge for the magnetic field 
B (c being the velocity of light and m the mass of the 
electrons) , 



Hn = 



2m 



-Bx, 



(1) 



The solution to this Landau level problem is of course 
familiar. The energy eigenvalues and eigenfunctions of 
an independent electron are 



en,k — fit^c 



1pn,k 



hck\ 



(2) 



where the frequency uic — — and the degeneracy of each 
energy level is 



= 2- 



(3) 



where the total flux threading the system is <i> = BL^Ly 
and the flux quantum is <i>o — hc/e; the factor of 2 is for 
spin. 

Even though momenta are not good quantum num- 
bers, it is still useful to depict the spectra on the two- 
dimensional kx — A:y-plane. The degenerate spectra, of 
degeneracy g, lie on concentric circles in this plane, sep- 
arated by hujc- Since momentum is no longer a good 
quantum number, the states are not located on specific 
points on the circle but can be viewed as rotating with 
frequency lOc- The total number of available states in a 
volume in the momentum space is unchanged in spite of 
the Landau quantization. In particular, if we denote AA 
to be the area between the concentric circles, then 



(2^)2 



(4) 



Imagine that the Fermi level, e^, at T = is situated 
on one such concentric level such that all states with 



energy E < ep are completely filled and all the levels for 
E > ep are completely empty. Then the total number of 
occupied states is eaxctly the same as the system without 
the magnetic field and the total energy per electron is also 
exactly the same. The area enclosed by the Fermi level, 
Alep), follows trivially from Eq. Q and is 



' (2^)^ 



(5) 



where N is the total number of particles. That is none 
other than the Luttinger sum rule. It is important to 
note that even though a magnetic field is never a small 
perturbation, the Luttinger sum rule is unchanged. 

The magnetic field corresponding to the ground state 
of a system with an integer number, n{ep), of Landau 
levels completely filled and all the rest completely empty 
will satisfy 



1 

Bn 



n{ep) 



27re 1 
he A{ep)' 



(6) 



Note that this ground state is an isolated nondegenerate 
state separated by a gap hcoc from the excited state. As 
we increase B, the quantized orbits are drawn out of 
the Fermi level, and sequentially pass through essentially 
identical set of nondgenerate isolated ground states. This 
of course results in the periodicity in the properties of the 
electron gas; the correction in the limit that n(ep) ^ 1 is 
negligible. Periodicity of course does not imply sinusoidal 
wave form and can contain higher harmonics. 

Now, fix Bn to a completely filled Landau level and 
turn on the electron-electron interaction. To all orders 
in perturbation theory, a nondegenerate isolated ground 
state will remain nondegenerate and therefore the se- 
quence of states corresponding to fully filled Landau lev- 
els as a function of the magnetic field will be the same, 
as in the noninteracting case. The periodicity is there- 
fore unchanged and is determined by the enclosed area 
A{ep), which in turn is fixed by the Luttinger sum rule. 

The argument will clearly break down if electron- 
electron interaction drives a quantum phase transition 
as a function of the ineteraction strength. A greater like- 
lihood for this happening is when the Landau level is 
partially filled and degenerate. Indeed, even for higher 
Landau levels, we know that a zoo of density wave states 
is a possibility,'^^ but Kohn's theorem should be robust 
for fully filled Landau levels. Note that Kohn's theo- 
rem makes no statement about the amplitude of oscilla- 
tions and is also silent about the waveform of the oscilla- 
tions. In general the periodicity, when Fourier analyzed 
will contain harmonics, and the harmonic content of the 
nonrelativistic case will be different from the relativistic 
case. 



B. Crystalline system: Dirac fermions in graphene 

We now extend Kohn's theorem to graphene. It is 
known from experiments that SdH oscillation frequen- 
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cies are in excellent agreement with the noninteracting 
system, as though electron-electron interaction plays no 
role. This is as it should be if the theorem holds. There 
are two basic issues that need to be dealt with: the crys- 
talline potential that breaks translational invariance and 
disorder in graphene samples. In this subsection we will 
discuss the former and will leave the latter to the follow- 
ing section. 

As before, first consider the non-interacting problem. 
The low energy spectrum for which the quasiparticles are 
well described by continuum Lorentz invariant Hamilto- 
nians based on two inequivalent nodes K = (|^, ) 

and K' = ( |^ , — ) in the tight binding description 
of the graphene lattice, with a the lattice spacing (cf. 
below) Near the Dirac point at K the energy eigen- 
functions are given by 



(7) 



where cr are the Pauli matrices and wf is the Fermi ve- 
locity. In the plane wave basis, the fc-space Hamiltonian 
Hk = hvpc-ii. Similarly, for K' , Hk' — hvpcr* -k. The 
spectra for each energy is two-fold degenerate, ignoring 
spin. 

The exact energy eigenvalues for this two-component 
Weyl fermion problem are of course trivial, and the corre- 
sponding eigenvalues and eigenfunctions in the presence 
of a perpendicular magnetic field are well known. While 
the energy eigenvalues differ from the corresponding non- 
relativistic Landau level problem, the spinor eigenfunc- 
tions are easily constructed from the corresponding non- 
relativistic problem. However, we will eschew the ex- 
act solutions and only take advantage of the validity of 
the continuum Hamiltonian to formulate the semiclassi- 
cal dynamics for large Landau levels that leads to the 
correct quantum oscillation frequencies. This is exactly 
what is necessary to demonstrate Kohn's theorem for the 
general case of electrons in a crystalline environment. We 
want to show that the Onsager quantization condition 
follows from the semiclassical quantization of the Landau 
orbits as the long wavelength low energy Hamiltonian is 
a valid description, regardless of whether not the descrip- 
tion is in terms of the relativistic Weyl Hamiltonian, as 
for graphene. 

In general, for the relativistic quasiparticles the semi- 
classical equations of motion of a Bloch electron are dif- 
ferent due to the presence of a nontrivial Berry curva- 
ture.'i^ These are 

r = ^VkS„(k)-kxJI„(k), (8) 



/ik = -eE - -r X B, 



(9) 



where E and B are the electric and magnetic fields re- 
spectively. The Berry curvature is defined as 



where 



n^{\i) = Vk X ^„(k), 



yl„(k) = i(u„(k)|VkUa(k)) 



(10) 



(11) 



is the Berry vector potential. The periodic part of the 
Bloch wavefunction is denoted by ^^(k), where a is the 
band index. 

When E = and B = Bz, following Chang and Niu,'^^ 
it can be shown that these set of equations of motion 
lead to the the following semiclassical quantization rule 
for the area of the orbit in momentum space 



Ak^^<l>ikx dk) 



2-KeB 

he 



(n-t-7). 



and the Maslov index 



r 

2^' 



with 



r = <j> Aa{k) ■ dk. 



(12) 



(13) 



(14) 



the Berry phase for the orbit. 

This phase easy is to compute for graphene because 
the eigenspinors for K and K' can be chosen to be 

^±,K - ^ ( ) , U±,K' = ^ ( ^gi^, ) (15) 

where ip\^ ~ tan^ ^ {ky / k^) . A simple computation shows 
that 



AK{k) = i{u±^K{k)\V^u±,K{k)) = - — , (16) 

AK'{k) = z(M±,x'(k)|VkU±.K'(k)) - ^. (17) 

where 9 is the azimuthal angle in the kx — ky plane and 
k = y/fc^ -I- ky. The Berry phase F is therefore ±7r. 
The semiclassical quantization formula is then 

Ak = -^—n. 18 
nc 

This implies once again that the magnetic field corre- 
sponding to a fully filled Landau level at the Fermi energy 
is given exactly by Eq. ([6]). The only remaining subtlety 
now is to note that the area of the Fermi pocket is 



nki 



2TTeB 
he 



(19) 



from which it immediately follows that the Fermi energy 
for the relativistic Dirac spectrum is given by 



ep = hvpkF = y^\/2n(eF). 



The magnetic length 



^hc/{eB). 



(20) 



(21) 



The appearance of this Berry phase F — ±tt for massless 
nodal quasiparticles with linear spectru m is a well known 
topological property of the Hamiltonian.'i^Iill 
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From here on the argument proceeds identically to the 
non-relativistic case, the only difference being the spac- 
ing between the Landau level, which is not constant as 
a function of n. As we turn on the electron-electron in- 
teraction, the energy levels are drawn out of the Fermi 
level in exactly the same sequence as the non-interacting 
case. Therefore the frequency of dHvA is unchanged. To 
the extent that the low energy continuum theory ade- 
quately describes graphene, there should be no effect of 
electron-electron interaction on the oscillation frequency. 

Due to the relativistic spectrum of massless quasiparti- 
cles in (2-1- l)-dimensions the density of states vanish lin- 
early at the nodal points, and the Coulomb interaction is 
not screened. This leads to strong long range interaction 
between the quasiparticles in graphene. In the presence 
of a moderate nonzero chemical potential the interaction 
will be screened, however. For small chemical potentials 
long range interactions can lead to broken symmetries 
and hence to a failure of Kohn's theorem. 

The nodal relativistic spectra can also arise from a 
condensate. An example is nodal fermionic quasiparti- 
cles of a particle-hole condensate in I = 2 angular mo- 
mentum channel, as in a singlet d-density wave (DDW), 
staggered flux phase, or an orbital antiferromagnet.^® 
This linearized Dirac fermion theory is valid for momen- 
tum small compared to inverse lattice spacing and en- 
ergy small compared to the bandwidth. For a wide range 
of the chemical potential, small compared to the band- 
width, one can still use the linearized continuum theory 
and results identical to those above hold. 



III. EFFECTS OF DISORDER 



A. Landau levels in graphene 



M-1 




N-1 N 



FIG. 1: (Color online) Graphene lattice. Lattice sites belong- 
ing to the same shaded horizontal or vertical stripes share 
the same indices i or j. The solid circles correspond to the 
A sublattice and the open circles to the B sublattice. The 
three nearest neighbor vectors joining the two sublattices are 
5k, fc = 1,2,3. 



eycomb lattice of dimension x Lj, in a perpendicular 
magnetic field, as shown in Fig. [T| which is 



n 

3 



(22) 



n k=l 



where the summation n ranges over all unit cells, and 
Cn,A, Cn,B are the fermionic annihilation operators in 
the unit cell n for the sublattices A and B, respec- 
tively. The spin degrees of freedom are omitted, as we 
assume that the magnetic field is sufficiently strong to 
completely polarize them. In principle, disorder can take 
various formSj^^^ but here we shall consider the on-site 
energy e^ A and Cn.s to be random variables uniformly 
distributed in the range [—gv/2,gv/2] corresponding to 
potential disorder. We choose the hopping matrix ele- 
ment tn,k = t = 1.0, providing a natural energy scale. 
The phases og^fc are chosen such that the magnetic flux 
per hexagonal plaquette, 4>, is l/Q, in units of the 
flux quanta (j)o = hc/e. We choose a gauge such that 
o-n.i ~ T^i/Q for the vertical bonds in slice i as in Fig. [T| 
and On, 2 = an,3 — 0. 

The Hall conductance Cxy can be computed by impos- 
ing periodic boundary conditions in both directions of the 
system and diagonalizing the Hamiltonian ( 22 1 to obtain 



a set of energy eigenvalues Ea and eigenvectors \a) for 



1, 



,L.r. X Ly. Then from the Kubo formula 



[201 



E 



{a\vy\/3){/3\v^\a) - {x <^ y) 
(23) 

where — [H,x\/ih is the velocity operator along the 
X direction and similarly for Vy. Note that the bonds in 
Fig. [l] that are not parallel to y direction contribute to 
both Vx and Vy. The summation corresponds to all states 
below and above the energy E. Finally, the expression is 
disorder averaged. 

The longitudinal conductance ctxx is studied using the 
well developed transfer matrix method. Consider a quasi- 
ID system, ^ Ly = 2M with a periodic boundary 
condition only along the y direction, where M denotes 
the number of unit cells in a slice along y direction. Let 
^'i (V'l,!, ^i.2, ■ ■ • , V'i,2M)"^ be the amplitudes on the 
slice i for an eigenstate with a given energy E; then the 
amplitudes on the successive slices are related by the ma- 
trix multiplication: 



Vr\E~H,) 



1 

^ ' '(24) 
diagonal matrix with elements 
representing the hopping matrix 
elements connecting the slices i and i + 1, and Hi is 
the Hamiltonian within the slice. All postive Lyapunov 
exponents of the transfer matrix^^ > 72 > . . . > 72Afi 



where Vi 



IS a 

, ti,2M) 



Consider the tight-binding Hamiltonian, iJ, on a hon- are computed by iterating Eq. (24 1 and performing 
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orthonormalization regularly. The convergence of this 
algorithm is guaranteed by the well known Osledec 
theorem J22l The cond uctance pe r square, cTxx, is given by 
the Landauer formula^^2l2ll2ll2ll(note the factor of \/3 in 
the argument of cosh due to the honeycomb lattice): 



2 2Af -, 

^ it cosh2(2V3M7,)' 



(25) 



The localization length in the quasi-lD system with Ly 
2M is given by Am = 1/72M- 



-9-8-7-6-5-4-3-2 -1 n=0 1 2 3 4 5 6 7 8 9 
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hopping parameter, t, which is set to unity. Each level is 
broadened due to disorder and finite size effects. Because 
the spacings between successive Landau bands become 
smaller as \n\ increases, the overlap between the neigh- 
boring bands increases, resulting in the overall parabolic 
shape of the background (Txx- The Hall conducivity a^y 
jumps by 2e^//i every time the energy passes through a 
Landau level, corresponding to the two fold valley degen- 
eracy in this system. 




FIG. 2: (Color online) The Hall conductance (Jxy (plateaus) 
and the longitudinal conductance CTxx as a function of en- 
ergy E. The magnetic flux through the hexagonal plaquette, 
(j) — 1/2000 in units of the flux quantum (j>o = hc/e and 
the potential disorder gv = 0.5,0.05,0.01,0.001 from top to 
bottom. For a^x, M was chosen to be 48, and, for a^y, the 
system size was chosen to be 4000 x 10 for the exact diago- 
nalization. In physical units the magnetic fleld corresponds 
to B ~ 40T. The vertical dashed lines indicate the locations 
of the Landau levels at En = sgn(n)(\/37r|n|/1000)^''^ for the 
parameters considered here. There are virtually no shifts of 
the positions of the Landau levels except for gv = 0.5 for 
which it is minimal. 

The results are shown in Fig. |2] where the Hall 
conductance a^y and the longitudinal conductance ctxx 
are computed as a function of energy E with param- 
eters (j) = 1/2000 and for potential disorder gy — 
0.5,0.05,0.01,0.001. The longitudinal conductance (Txx 
peaks almost exactly at the Landau levels 



E„ = sgn{n)h ( V2 



-VF 



where 



v-p = 3ta/2h 



(26) 



(27) 



is the fermi velocity, a being the lattice spacing. The 
only exception is the case gv = 0.5, as shown in Fig. |3] 
for which there is a minimal shift. Recall that this is 
very large disorder as the energy scale is in terms of the 



FIG. 3: (Color online) The location of the peaks in Fig. [5] 
Note that the slope is essentially the same as the noninter- 
acting system except for gv = 0.5 for which it is slightly 
different. 



B. SdH Oscillations: analytical results in the 
self-consistent Born approximation 

Shubnikov-de Haas oscillations in graphene can be cal- 
culated numerically by the method described in the previ- 
ous subsection. However, for realistic fields, it is difficult 
to control the accuracy because of the essential singu- 
larity corresponding to the magnetic field corresponding 
to all quantum oscillation phenomena. Moreover, the 
dependence on the physical parameters are not particu- 
larly transparent. For this reason we adopt an analytical 
SCBA developed by Ando.E^l 

In the presence of disorder, the imaginary part of the 
self energy at the Fermi level is non-zero. Hence, Lut- 
tinger's theorem and Kohn's argument are not immedi- 
ately applicable. However, for weak disorder, the quasi- 
particle lifetime at the fermi energy, t{£f) = r (hence- 
forth by T we will mean T(ep)), can be long and the 
Fermi surface can be reasonably well defined within the 
uncertainty h/r. If the Fermi energy ep is very large, 
we can use h/epT as a small parameter. Due to disor- 
der the Landau levels are broadened into bands and in 
the limit of overlapping Landau levels (wcT < 1) explicit 
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calculations using SCBA show that the period of oscilla- 
tions is unchanged to an excellent approximation. The 
correction terms to the period is of the order \hl{epT)Y. 
For unitary scatterers, the shift is much larger, - but it is 
doubtful that such strong potential scattering is relevant 
to grapheneP 

Since the gap between the relativistic Landau levels 
decreases for higher Landau levels, disorder will have 
stronger effect in contrast to the equally spaced non- 
relativistic Landau levels. Hence, in the semiclassi- 
cal limit of higher Landau levels, the damping will be 
stronger for the relativistic case. For the sake of sim- 
plicity consid er th e short range impurity potential. Fol- 
lowing Ando,'^2l23 I; Jig self-consistency equation for the 
self-energy E(e) in SCBA can be written as 



S(6) 



Atteft 



E 



1 



e-En- S(e) ■ 



(28) 



Note that we have introduced a cutoff Nc, which is of 
the order of the bandwidth. It determines the limit of 
the applicability of the linearized Dirac spectra. Let 

e-S(e) - X{e)^X' + iX"{e), (29) 
a(e) - x'ie)/ihio), (30) 
6(e) = X"(e)/(M, (31) 

where iv = ^/2vf/Ib ■ Ignoring the arguments of a and 
h, the self-consistency equation becomes 



where « t&n~^{l%/2epT), and y = 1 — 3{h/2epT)'^ . 
We have defined a mass parameter m* — ep/vp, and 
Lu* — eB/m*c, the cyclotron frequency of a hypothetical 
nonrclativistic system. Here the frequency F is 



and the area of the Fermi pocket is given by 



(35) 



(36) 



The condition a '3> b implies (fi./2e^r) <C 1, which justi- 
fies the weak disorder approximation. When {h/2epT) ~ 
1, the SCBA is insufficient to treat the randomness cor- 
rectly. Hence, in our calculation we assume {h/2epT) <C 
1. 

Using the Kubo formula the longitudinal dc conduc- 
tivity axx becomes 



°° df(e) 
oe 



(37) 



where /(e) is the Fermi-Dirac distribution function and 
the kernel K^x i 



HlOc 



a — ib 



2TrepT 



E 



(a + ib) 



1 1 



(a -|- ib)^ — m 2 a + ib 



(32) 

The sum on the right hand side can be calculated using 
the Poison summation formula described in Appendix [A] 
and, to leading orders in {h/epr), a and b are given by 



e 

huj 



1 



h 

2epT 



— e sin ( — — 



epT 
e 

2epLUT 



B 



l + 2e~ 



T^^^p' cos ( — |— 



(33) 
(34) 



Kxxie) 



J2^9mie + i0)^gm+iie + i0) (38) 



m=0 



In the above equation 



5m (e) 



e- S 



(e-E)2-m(M^ 



(39) 



At r = 0, axx{T = 0) = Kxxiep)- After using the Pois- 
son's summation formula the conductivity can be written 



as 



axx ~ cr 



l-,2e-'^/--|cosr^ 



(c^,*r)2-l / 2^^^(e;^) 
cos y- 



(c.*r)2 + 1 



B 



3h . / 27TF(ep) 

- sm y- 



2epT 



B 



(40) 



where 



(^;r)2 + l' 



and CTq = An^e^T/Trm* ; is the density of quasiparti- 
cles corresponding to one of the valleys for a given spin 
direction. 

(41) 

^ ' In contrast to the non-relativistic Landau levels, the 
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frequency of the oscillation is changed by a factor y and 
there is a small phase shift in the oscillations. Hence, 
the robustness of the oscillation frequency and the phase 
of the oscillations do depend on on disorder but only 
very weakly. In the SdH experiments, this phase shift 
has not been observed and this can be attributed to the 
weakness of the disorder. When the disorder is weak 
i.e, {h/2eFT) <C 1, y « 1, </> 0, and the conductivity 
becomes 



1 



-tt/o. 



1 



■ cos 



27rF(eF) 
B 



(42) 



If a-'*r < 1 the oscillations will be heavily damped. 

Although this expression appea rs ide ntical to the cor- 
responding nonrelativistic formula,'^^^ it is actually dif- 
ferent. The definition of to* depends on the m* that we 
have defined, which in turn depends on ep- Therefore, 
the oscillations will be damped more strongly for higher 
Landau levels. The reason behind this is the smaller 
gaps between the higher Landau levels, as compared to 
the nonrelativistic case. 



IV. CONCLUSIONS 

In conclusion, we reiterate that the positions of the 
Landau levels, as well the frequencies of SDH oscillations, 
are remarkably robust with respect both impurity scat- 
tering and electron-electron interaction (as long as the 
chemical potential is not too close to the Dirac points) . 
This seems to be consistent with experiments .'21 In the 
presence of both interaction and randomness the Lut- 
tinger's many body formalisnP becomes applicable, even 
for two dimensions in the presence of disorder and/or 
thermal broadening. From the above considerations, if 
the disorder is small enough {h/tpT ^ 1), the oscilla- 
tion frequencies can be expected to be proportional to 
the true Fermi surface area of the pure interacting prob- 
lem. In principle one can attempt to treat the interac- 
tion within a Hartree-Fock (HF) approximation to esti- 
mate the effects of the interaction on the amplitude of 
the oscillations. But, there is no reason to trust the HF 
results. The correlation energy contributions can have 
equally important effects on the amplitude, in particu- 
lar, on the effective mass parameter. One also needs to 
account for the disorder induced vertex corrections for 
the self-energy contributions from the interaction. Such 
a calculation will be tedious but important to understand 
the effects of inelastic scattering rates and we shall rele- 
gate such detailed calculations for a future publication. 
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APPENDIX A: SCBA 



Consider the real and the imaginary parts of the Eq. 
(I32| : 



47rei?T / e 

I a 



1 



= 2aS 
= 2S 



+ 1 



(Al) 
(A2) 



where 



'5- = E 



+ &^ ± TO 



pN, 



Jo 



■Jj (a^ - &2 _ „i)2 _|_ 40,2^2 

oo 

dm[l H- 2 cos(27rA:m)] 



(A3) 



k = l 



(a2-62-TO)2-h4a262 
Now choosing m ~ + b'^ ^ z, we obtain 



(A4) 



5+ = 



dz[l + 2 V cos(27rfc(a2 - b^ + z))] 

k=l 

.2 



2a2 + z 



(A5) 



z2 -I- 40^62 ' 

dz[l + 2 ^ cos(27rfc(a2 - b^ + z))] 

k=l 



' z2 + 4a2^2 ■ 



(A6) 
(A7) 



Because Nc ^ a(e)^ ^ K^)'^ f™' ^ ^ ^f, we can let the 
upper limit tend to oo and the lower limit to —oo. Note 
that after this change of limits the terms that are odd in 
z vanish. We find the following analytic expressions: 



h h ^ 



fe=l 



X COS 



27rfc(a^-62) + tan~^ 



(A8) 



5_ = !^ + ?JIv/a^£e-4'^'="^ 

k=\ 



a a 



X sm 



27rfc(a2 - 5^) +t8 



(A9) 



Substituting S± in the Eqs. ( Al land ( A2 ), retaining only 



the first harmonic and onl y te rms to l eading order in 
{h/eprY we arrive at Eqs. (33l and (34l. 
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The conductivity kernel K^^ in Eq. ( 38 ) can be ex- 
pressed as 



No 

E 



a +b'^ + m 



■K^h ^ id?- - IP- - niY + 40^62 

m— 

+ 6^ + m + 1 



(AlO) 



{a^ - IP- -m- 1)2 + 4a262 
Now using Poisson's summation formula we get 

— -[/i + y cos(2^A:(a2 _ \?))h{k) 

k=l 

oo 

- sin(27rfc(a2 _ b^))l3{k)] (All) 



K — 



where 

^2(fc) = 



fe=l 



°° (2a^ + z)(2o^ + z + l) 

^[z2+4a252][(^+ 1)2 + 4^2^2] 



(A12) 



cos(27rfcz)(2a^ + z)(2a^ + ^ + 1) 
, [z2+4a2fe2][(z+ 1)2 +4^252] ^ ) 



We can clearly see 2/i = /2(A: = 0), and all the integrals 
can be evaluated by performing a single integral 



(2a2 + z)(2a^ + z+l) 
[z2+4a262]p + 1)2 +4^252] 



(A15) 



The integrand has simple poles ai z — ±2iab and z — 
— 1 ± 2iab, and we evaluate by closing the contour in 
the upper half-plane. It turns out that = Iml^ = 
and 



hik) = Reh 



-47Tkab 



1^ + 9 
1 + 160252 



(A16) 



Now sub stituting the values of these integrals in 
Eq. (All I, and also retaining only the first harmonic of 



the oscillations we obtain 



K. 



8e^ab + &2 



1-H2e" 



4:7Tab 



TT^h 1 + 16a2&2 ^ 

X cos(27r(a^ - 6^)) 



(A17) 



^ sin(2^fcz)(2a^-Hz)(2a^H-z-H) 
^ y_ [z2+4a2fe2][(,+ i)2+4a262] (^14) 



Now substituting the expressions for a and b and keeping 
terms upto 0{h/eF'rY we obtain Eq. (40 1. 



^ D. Shoenberg, Magnetic oscillations m metals (Cambridge 

University Press, Cambridge, 1984). 
^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 

A. A. Firsov, Nature 438, 197 (2005). 
^ N. Doiron-Leyraud, C. Proust, D. LeBoeuf, J. Levallois, 

J.-B. Bonnemaison, R. Liang, D. A. Bonn, W. N. Hardy, 

and L. Taillefer, Nature 447, 565 (2007). 

* For a recent perspective, see S. Chakravarty, Science 319, 
735 (2008). 

^ W. Kohn, Physical Review 123, 1242 (1961). 

® J. M. Luttinger, Physical Review 121, 1251 (1961). 

S. Curnoe and P. C. E. Stamp, Physical Review Letters 

80, 3312 (1998). 

* J. M. Luttinger, Physical Review 119, 1153 (1960). 

^ S. G. Sharapov, V. P. Gusynin, and H. Beck, Physical 

Review B 69, 075104 (2004). 
^° N. M. R. Peres, F. Guinea, and A. H. Castro Neto, Physical 

Review B 73, 125411 (2006). 
" L M. Lifshits and A. Kosevich, Zh. Eksp. Teor. Fiz. 29, 

730 (1955), [Sov. Phys. JETP 2, 636 (1956)]. 

T. Ando, .Journal of the Physical Society of Japan 37, 1233 

(1974). 

A. A. Koulakov, M. M. Fogler, and B. L Shklovsku, Phys- 
ical Review Letters 76, 499 (1996). 

A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, and A. K. Geim, The electronic properties of 
graphene (2008), http://arxiv.org/abs/0709.1163. 
M.-C. Chang and Q. Niu, Physical Review B 53, 7010 



(1996). 

G. E. Volovik, The universe in a Helium droplet (Cam- 
bridge University Press, Cambridge, 2003). 
J. L. Maries, F. Guinea, Maria A. H. Vozmediano, Phys. 
Rev. B, 75, 155424 (2007). 

S. Chakravarty, R. B. Laughhn, D. K. Morr, and C. Nayak, 

Phys. Rev. B 63, 094503 (2001). 
^® X. Jia, P. Goswami, and S. Chakravarty, Physical Review 

Letters 101, 036805 (2008). 
^° D. J. Thouless, M. Kohmoto, M. P. Nightingale, and 

M. den Nijs, Phys. Rev. Lett. 49, 405 (1982). 

B. Kramer and M. Schreiber, in Computational Physics, 

edited by K. H. Hoffmann and M. Schreiber (Springer, 

Berlin, 1996), p. 166. 

V. Oseledec, Trans. Moscow Math. Soc. 19, 197 (1968). 
23 D. S. Fisher and P. A. Lee, Phys. Rev. B 23, 6851 (1981). 
2" H. U. Baranger and A. D. Stone, Phys. Rev. B 40, 8169 

(1989). 

2^ B. Kramer and A. MacKinnon, Rep. Prog. Phys 56, 1469 
(1993). 

2^ D. N. Sheng and Z. Y. Weng, Europhys. Lett. 50, 776 
(2000). 

Y. Zheng and T. Ando, Physical Review B 65, 245420 
(2002). 



